Quantum quenches and thermalization in one-dimensional fermionic systems 
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We study the dynamics and thermalization of strongly correlated fermions in finite one-dimensional lattices 
after a quantum quench. Our calculations are performed using exact diagonalization. We focus on one- and 
two-body observables such as the momentum distribution function [n(k)] and the density-density structure 
factor [N(k)], respectively, and study the effects of approaching an integrable point. We show that while the 
relaxation dynamics and thermalization of N(k) for fermions is very similar to the one of hardcore bosons, 
the behavior of n(k) is distinctively different. The latter observable exhibits a slower relaxation dynamics in 
fermionic systems. We identify the origin of this behavior, which is related to the off-diagonal matrix elements 
of n(k) in the basis of the eigenstates of the Hamiltonian. More generally, we find that thermalization occurs 
far away from integrability and that it breaks down as one approaches the integrable point. 

PACS numbers: 03.75.Ss, 05.30.Fk, 02.30.1k, 67.85.Lm 



I. INTRODUCTION 

In recent years, the study of the nonequilibrium dynamics 
of isolated quantum systems has attracted a great deal of at- 
tention. The main motivation for these studies lies behind 
the spectacular success that experimentalists have achieved 
in trapping and manipulating ultracold quantum gases. This 
has allowed them to, for example, load ultracold bosons in 
optical lattices and study the collapse and revival of the mat- 
ter wave interference after quenching the interaction strength 
from deep in the superfluid regime into deep in the Mott in- 
sulating regime [1], and to observe the damping of the dipole 
oscillations when the center of mass of the gas was displaced 
away from the center of the trap S IH . 

Since these atomic gases are trapped, cooled, and manip- 
ulated in a very high vacuum, i.e., in no contact with any 
kind of thermal or particle reservoir, a fundamental question 
that arises is whether after a sudden perturbation the system 
will relax to thermal equilibrium (an equilibrium in which ob- 
servables are described by standard statistical mechanical en- 
sembles). The usual assumption is that thermalization occurs 
in general. However, recent experiments in one-dimensional 
(ID) geometries (created by a deep two-dimensional optical 
lattice) have failed to observe relaxation to a thermal distribu- 
tion after a quench |5[]. The absence of thermalization can be 
understood in the very special regime in which the system is 
at an integrable point [6]. Interestingly, in Ref. [5] the gas was 
perturbed away from integrability and thermalization did not 
occur. Other experiments in ID geometries (generated by an 
atom chip) have reported the indirect observation of thermal- 
ization [7]. In the latter case, the transverse confinement was 
not as strong as the one induced by the optical lattice in Ref. 
|Ht]. The question of whether strict one dimensionality could 
affect the outcome of the relaxation dynamics has not been 
fully addressed experimentally. 

The particular case in which the system is at an inte- 
grable point lends itself to combined analytical and numer- 
ical studies, which have allowed theorists to show that in 
that case the relaxation dynamics after a quench lead to non- 
thermal equilibrium distributions of few-body observables 
|@]. In addition, a generalization of the Gibbs ensemble, 



which takes into account the conserved quantities that make 
the system integrable, has been shown to successfully de- 
scribe experimental observables after relaxation [6, 8]. Fur- 
ther analytical and numerical works have analyzed the rele- 
vance of the generalized Gibbs ensemble to different mod- 
els and observables and addressed its limits of applicability 

aiiaimEiiaHEEiiziiiia]. 

When the systems are far away from integrability, for exam- 
ple in higher dimensions, the expectation is that they should 
thermalize. Recent numerical work has shown this to be the 
case for bosons in two dimensions, studied with exact diago- 
nalization B19I1 . and for fermions studied within the dynamical 
mean-field theory (DMFT) approximation lEoll . In addition, 
thermalization in the isolated two-dimensional bosonic case 
could be understood within the eigenstate thermalization hy- 
pothesis (ETH) proposed by Deutsch 02111 and Srednicki 12211 . 
in which the individual eigenstates of the many-body Hamil- 
tonian exhibit thermal behavior 1U9I1 . 

In one-dimensional nonintegrable systems, the situation is 
a bit more subtle as many models can be tuned (by chang- 
ing some Hamiltonian parameters) to be arbitrarily close to 
integrable points. An early study of the relaxation dynamics 
in fermionic systems when breaking integrability, using time- 
dependent density renormalization group (tDMRG) 12311 . con- 
cluded that thermalization does not occur even if the system is 
perturbed away from the integrable point. Another early work 
in ID, also using tDMRG, studied the relaxation dynamics of 
the (nonintegrable) Bose-Hubbard model when quenching the 
system across the superfluid-to-Mott-insulator transition Il24ll . 
In that case the authors concluded that thermalization occurs 
in some regimes but not in others. 

More recently, we have performed a systematic study of 
thermalization in finite ID systems of hardcore bosons. Our 
results indicate that far away from integrability the system 
does thermalize [|25ll . However, thermalization breaks down 
as one approaches the integrable point. An important ques- 
tion that still needs to be answered is what happens with the 
point at which thermalization breaks down as one approaches 
the thermodynamic limit. It may either move toward the in- 
tegrable point or toward some point away from integrability. 
More powerful numerical techniques or analytical approaches 
may be required to answer that question, which will not be 
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addressed here. 

In this work, we extend the analysis in Ref. B25I1 to the 
fermionic case, which allows us to discuss some of the is- 
sues that remained open in Ref. ll23ll . In particular, the time 
scales required for different observables to relax to a station- 
ary "equilibrium" distribution, and the description of the ob- 
servables after relaxation. We also study how the closeness to 
an integrable point affects the dynamics and thermalization af- 
ter a quantum quench. As in earlier works lfl9ll25ll . we show 
that the breakdown of thermalization for our observables is 
directly linked to the breakdown of ETH. 

The presentation is organized as follows. In Sec. [II] we 
introduce the lattice model to describe ID fermions with 
nearest- and next-nearest-neighbor hopping and repulsive in- 
teractions. We define our observables of interest and briefly 
discuss our numerical approach. The nonequilibrium dynam- 
ics of this model is presented in Sec. [Ill] Section[lV]is devoted 
to the analysis of statistical mechanical approaches to describe 
observables after relaxation. In Sec.[V] we justify these statis- 
tical mechanical approaches, or their failure, within ETH. An 
explanation of why the momentum distribution function of the 
fermions has a distinctively slow relaxation dynamics is pre- 
sented in Sec.[Vl]in terms of the off-diagonal elements of this 
observable in the basis of the eigenstates of the Hamiltonian. 
Finally, the summary and outlook are presented in Sec. IVHI 

H. HAMILTONIAN AND OBSERVABLES 

The Hamiltonian of spin-polarized fermions in a ID lattice 
with periodic boundary conditions can be written as 

^ = E {"* +H.c.) +v[n i - (n i+1 - 

- t' [flf i+2 + H.c.) + V Li - (n i+2 - | . (1) 

where the fermionic creation and annihilation operators at site 
i are denoted by f! and respectively, and the local-density 
operator by hi — fjf^. In Eq. (fl~|i, the nearest- and next- 
nearest-neighbor hopping parameters are denoted by t and t', 
respectively, and the nearest- and next-nearest-neighbor inter- 
actions are denoted by V and V', respectively. In our study, 
we only consider repulsive interactions (V, V > 0), and the 
number of lattice sites is denoted by L. 

In Ref. rt25ll . we have already studied a very similar Hamil- 
tonian for hardcore bosons (bosons with an infinite on-site re- 
pulsion), in which case 

H b = J2 j-t (btb l+ i + H.c.) + V [n\ - 1) [n\ +l 

- H (btU+2 + H.c.) + V (n\ - [n\ +2 - } . (2) 

where the hardcore boson creation and annihilation operators 
at site i are denoted by b\ and b it respectively, and the local- 
density operator by h\ — SjS^ The parameters t, V, t' , and 



V' in Eq. (O have exactly the same meaning for hardcore 
bosons as for the fermions in Eq. (Q~|i. For hardcore bosons, 
the creation and annihilation operators commute as usual for 
bosons 

M]] = Mi] = ®,&5] = > for ^3, (3) 

but on the same site, the hardcore bosons operators satisfy 
anticommutation relations typical of fermions, 

{bJt} = l, bf = bi = Q. (4) 

These constraints avoid double or higher occupancy of the lat- 
tice sites. 

Similar to the hardcore boson case studied in Ref. ll25ll . the 
Hamiltonian ([U is integrable for t' = V' = 0. In this section, 
we restrict our analysis to systems with 1/3 filling, and study 
the effects of finite, but small, values of t' = V ^ when 
one departs from the integrable point. In this case, the ground 
state of our Hamiltonian is always a Luttinger liquid and no 
metal-insulator transition occurs in the system 1126I1 . 

For t' = 0, the fermionic Hamiltonian ([Til can be mapped 
onto the identical Hamiltonian for hardcore bosons (O (up to 
a boundary effect), i.e., the fermionic operators in Eq. ([T) just 
need to be replaced by operators describing hardcore bosons. 
This can be done using the Jordan- Wigner transformation 
1I27I1 . and implies that in the thermodynamic limit both sys- 
tems have the same spectrum and identical diagonal correla- 
tions. However, off-diagonal correlations and related observ- 
ables, such as the momentum distribution function, are very 
different for fermions and bosons. A finite value of t' still al- 
lows for a mapping of the fermionic Hamiltonian dTJ onto a 
Hamiltonian of hardcore bosons. However, the new hardcore 
boson Hamiltonian for t' ^ will be different to the one in 
Eq. (f2]l. The mapping introduces additional operators in the 
term proportional to t'. In addition, in the fermionic case, a 
finite value of t' breaks the particle-hole symmetry present for 
t' = and present for any value of t' in the hardcore boson 
case. 

The study of the nonequilibrium dynamics and thermody- 
namics of these systems is performed using full exact diago- 
nalization of the Hamiltonian dTJ. For our largest system sizes, 
eight fermions in 24 sites, the total Hilbert space has dimen- 
sion D = 735 471. We take advantage of the translational 
symmetry of the lattice, which allows us to block diagonalize 
the full Hamiltonian while the dimension of the largest block 
(in k space) is Dk = 30 667. 

In this work, we focus our analysis on two observables. The 
first one in the momentum distribution function 

n(*) = ^E e "* ( *" J ' ) ^/j' (5) 

id 

which is the Fourier transform of the one-body density matrix 
(f>ij = fjfj)- A s mentioned before, this observable behaves 
very differently for fermions as compared with the hardcore 
bosons to which they can be mapped. We should add that n(k) 
is usually measured in experiments with ultracold quantum 
gases via time-of-flight expansion. 
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The second observable of interest is the density-density 
structure factor 



N(k) = i^e-^Mj, 



(6) 



which is the Fourier transform of the density-density correla- 
tion matrix. Since we work at fixed number of fermions Nf, 
the expectation value of N(k = 0) is always (N(k — 0)) = 
iV 2 I L and, as usual, we set it to zero by subtracting that trivial 
constant. The structure factor can be measured in cold gases 
experiments by means of noise correlations |28il . 



HI. NONEQUILIBRIUM DYNAMICS 

We will restrict our study of the nonequilibrium dynamics 
to the case in which the system is taken out of equilibrium by 
means of a quench. Quenches are a special way to induce dy- 
namics by starting with an eigenstate of some initial Hamilto- 
nian H ini , and then instantaneously changing the Hamiltonian 
(at time r = 0) to some final time-independent Hamiltonian 

Hf in . 

As mentioned previously, our model is integrable for t' = 
V = 0. Since we are interested in studying the effect that ap- 
proaching an integrable point has on the dynamics and ther- 
malization of the system, we take the initial state to be an 
eigenstate of a system with t = t ini , V = Vi n i, t', V', 
and then quench the nearest-neighbor parameters t and V to 
t = tfi n , V = Vfm without changing t', V', i.e., we only 



change t mi ,V ini -> t fin ,V fin . This quench is performed 
for different values of t', V' as one approaches the integrable 
point t' = V = 0. 

Notice that our quenches do not break translational symme- 
try. For that reason, we always select the initial state to be one 
of the eigenstates of the initial Hamiltonian in the total k — 
sector. Since translational symmetry is preserved, only states 
with zero total momentum are required to calculate the exact 
time evolution of the system 



Mr)) 



,-iH fin r 



IVmi) =Y, C ^ lEaT \0 > C7) 



where \ip(r)) is the time-evolving state, \ipini) is the initial 
state, are the eigenstates of the final Hamiltonian with 
zero total momentum, energy E a , and C a — V'im}- The 
sum runs over all the eigenstates of the total k = sector. 
In this section, we perform the exact time evolution of up to 
eight fermions in lattices with up to 24 sites. This means that 
the largest total k = sector diagonalized had 30666 states. 

In contrast to classical systems, where one has to perform 
the time evolution in order to compute the long-time aver- 
age of any observable, in isolated quantum systems (where 
the time evolution is unitary) one can predict such time av- 
erages without the need of calculating the dynamics. Since 
the wave function of any initial state can be written in terms 
of the eigenstates of the final Hamiltonian \^ a ) as \ipini) = 
Col |^a}> one finds that the time evolution of the quantum 



expectation value of an observable O can be written as 

(0(r)) = <V(t)|6|V(t)> ^(ZCpeWo-^Ocp, 

a, 

where O a p are the matrix elements of O in the basis of the 
final Hamiltonian. This in turn implies that the infinite time 
average of the observable can be written as 



(o(r)> =o diag = ]T|cy 2 o QQ . 



(9) 



We have assumed that, as in the case of generic (noninte- 
grable) systems, the spectrum is nondegenerate and incom- 
mensurate. This means that if the expectation value of O re- 
laxes to some kind of equilibrium value (up to the recurrences 
that must occur if the system is isolated), that value should be 
the one predicted by Eq. (O. Following previous work, we 
prefer to think of this exact result as the prediction of a "diag- 
onal ensemble," where \ C a | 2 is the weight of each state within 
this ensemble 11191 12511 . 

One objection that may arise at this point is that even if 
the system is nondegenate and incommensurate, the spectrum 
may have arbitrarily close levels and it may take an unrealisti- 
cally long time for the prediction in Eq. Q to apply, i.e., that it 
may not be relevant to describe experiments. In order to study 
how the dynamics drives our observables of interest [n(k) and 
N(k)] toward the prediction of Eq. (0, we follow the scheme 
in Ref. [25]. We study the normalized area between n(k,r) 
and N(k,r), during the time evolution, and their infinite time 
average, i.e., at different times we compute 



and 



5n k (r) 



SN k (r) 



Efc K fc , T ) - n diag {k)\ 

Tldiag \k) 



E k \N(k,r)-N dtag (k)\ 



22k Ndiag 



(k) 



(10) 



(11) 



In Fig.Q] we show results for Sn k (T) and SNk(r) as a func- 
tion of time t for four different quenches as one approaches 
the integrable point. Besides undergoing the same change 
t^i = 0.5, Vini = 2.0 -> tfm = 1-0, Vfin = 1-0, these 
quenches have another very important property in common. 
The initial state for each was selected to be one eigenstate of 
the initial Hamiltonian in such a way that the effective tem- 
perature T of the system [29] is always the same (T = 2.0 
in Fig.Q]). Given the energy of the time-evolving state in the 
final Hamiltonian, which is conserved, 

E = (4>ini\Hfi n \lpini), (12) 

the effective temperature T ll29ll is defined by the expression 

E=^Tv[Hf ln e-"^' T Y (13) 

where 

Z = Tr|e-^ i " /T | (14) 
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is the partition function, and we have set the Boltzmann con- 
stant ks to unity, and tf in = 1 sets the energy scale in our 
system. 

In Figs. [If a) and [TJf), we compare the initial momentum 
distribution and structure factors with the predictions of Eq. 
© away from integrability (f = V = 0.32). (The re- 
sults for other values of t' , V are similar and not shown 
here.) The dynamics of those observables for four different 
values of t' , V' is presented in Figs.[Ttb)-02e) for 5nk(r) and 
in Figs. OlgHIlj) for <^Vfc(r). They show that the dynam- 
ics of n(k) and N(k) are different. Away from integrabil- 
ity (t' = V 1 ^ 0), SNk(r) quickly relaxes [in a time scale 
r ~ (tfin)^ 1 = 1.0] to a value ^0.02 and then fluctuates 
around that value. 8nk(r), on the other hand, slowly drifts to- 
ward 5nk{r) ~ 0.02-0.04 in a much longer time scale, orders 
of magnitude longer than the one seen for the relaxation of 
SN k {r). 
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FIG. 1: (Color online) Quantum quench tini = 0.5, Vi„i = 2.0 
-> t fin = 1.0, Vfm = 1.0, with t' ini = t' fin = t' and V( ni = 
Vj in = V' in a system with Nf = 8 and L = 24. The initial state 
was chosen in such a way that after the quench the system has an 
effective temperature T = 2.0 12911 in all cases. Given the energy 
of the initial state in the final Hamiltonian E = {il>%ni\H fi n \ipini), 
the effective temperature is computed following Eq. d 1 3b - (a) Initial 
and diagonal ensemble results for n(k) when t' = V' = 0.32. (b)- 
(e) Time evolution of 8n k for t' = V' = 0.32, 0.12, 0.04, and 0.0, 
respectively, (f) Initial and diagonal ensemble results for N(k) when 
t' = V = 0.32. (g)-(j) Time evolution of SN k for t' = V' = 0.32, 
0.12, 0.04, and 0.0, respectively. 

The behavior of SNk (r) in these spin-polarized systems is 
qualitatively (and quantitatively) very similar to the one ob- 
served for hardcore bosons in Ref. j25|]. The time evolution 
of dnk for the fermions, on the other hand, is different from 
the one seen for the hardcore bosons in Ref. ll25ll . For the 
latter systems, n(k) quickly relaxed toward the predictions of 
the diagonal ensemble, even at integrability. In Fig.[Hj), one 
can see that at integrability 8Nk(r) for the fermions also re- 
laxes quickly toward the predictions of the diagonal ensemble 
and then exhibits fluctuations between SNk(r) = 0.02 and 
SNk(r) — 0.06, while 5nk(r) in Fig. [He) exhibits very large 
oscillations between 5nk{r) — 0.02 and Snk(r) — 0.12. Ad- 
ditionally, we studied the dynamics of 5rik for a time scale ten 



times longer than the one depicted in Fig. [TJe) and observed 
exactly the same behavior. 

For hardcore bosons II25I]. we have shown that increasing 
the effective temperature 112911 (or the final energy E) de- 
creases the mean value of Snk(r) and 5Nk(r) after relaxation 
and also decreases the temporal fluctuations of those quanti- 
ties. In Fig. |2] we show results for quenches similar to the ones 
in Fig. [T] but with a higher effective temperature (T = 3.0). 
The comparison between Figs.|2]and[T]also reveals a reduction 
in the mean values of Snk(r) and SNk(r) and their fluctua- 
tions, which is particularly obvious at integrability. Still, dif- 
ferences remain between the relaxation dynamics of Snk(r) 
and 5Nk(r). Also, one can notice that at integrability the 
mean values of Snk(r) and SNk(r) and their fluctuations are 
in general larger than away from integrability, and they all de- 
crease the further away one moves from the integrable point. 
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FIG. 2: (Color online) Quantum quench tini ~ 0.5, Vini = 2.0 
-» tfi„ = 1.0, Vfin = 1.0, With t'ini = t'fin = l ' and V Li = 
V'fin = V in a system with Nf = 8 and L = 24. The initial 
state was chosen in such a way that after the quench the system has 
an effective temperature T = 3.0 ] 29] in all cases (see the caption 
of Fig. [TJ. (a) Initial and diagonal ensemble results for n(k) and 
t' = V = 0.32. (b)-(e) Time evolution of 5n k for t' = V = 0.32, 
0.12, 0.04, and 0.0, respectively, (f) Initial and diagonal ensemble 
results for N(k) andt' = V' = 0.32. (g)-(j) Time evolution of 6N k 
for t' = V' = 0.32, 0.12, 0.04, and 0.0, respectively. 

The improvement of the relaxation dynamics with increas- 
ing the effective temperature (energy) of the isolated state 
112911 . discussed above, can be related in part to the increase 
in the density of states in the final system with increasing en- 
ergy. This increases the number of eigenstates of the Hamil- 
tonian that participate in the dynamics making dephasing in 
Eq. ^ more effective, and hence reducing temporal fluctua- 
tions after relaxation. This can be better seen in Fig. [3] which 
clearly shows that the number of states with the largest val- 
ues of \C a \ 2 increases as the temperature decreases. Since 
J2 a \C a \ 2 = 1, this means that the lower the temperature the 
smaller the total number of states that participate in the dy- 
namics. This can be also seen in Fig. [3] where the number 
of states with \C a \ 2 > 10~ 5 and |C Q | 2 > 10~ 4 is larger for 
T = 3.0 [Fig.db)] than for T = 2.0 [Fig. Ha)] (notice the 
logarithmic scale in both axes). Overall, one can conclude that 
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dephasing becomes less effective with decreasing temperature 
and this leads to an enhancement of temporal fluctuations. 

A similar argument, based on the number of states partic- 
ipating in the dynamics, can help us understand why the av- 
erage values and temporal fluctuations of Snk (r) and 5N k (r) 
are in general larger after relaxation when one is at integra- 
bility, or close to it. Figure [3] also shows that for any given 
temperature, the total number of states that participate in the 
dynamics decreases as one approaches the integrable point. At 
integrability, this may be related to the presence of additional 
conserved quantities, which restrict the number of eigenstates 
of the final Hamiltonian that can have a significant overlap 
with the initial state. 
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FIG. 4: (Color online) Quantum quench tini = 0.5, Vini = 2.0 
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FIG. 3: (Color online) Number of states with \C a \ 2 greater than the 
value presented in the x axis, for an effective temperature (a) T = 
2.0 and (b) T = 3.0, and for the same quenches studied in Figs.Q] 
and[2] respectively. Here, Nf = 8 and L — 24. 

It would be desirable to study how the mean values and 
temporal fluctuations of Snk(r) and 5Nk(r) behave after re- 
laxation as one increases the system size. Unfortunately, the 
computational requirements of our approach increase expo- 
nentially with system size and a rigorous finite-size scaling is 
not possible. As a step in understanding how finite-size effects 
affect our results here, we show in Fig. HJresults for a smaller 
system with seven fermions in 21 lattice sites (with T — 3.0). 
By comparing those results with the ones in Fig. [2] one can 
see that, as expected, decreasing the system size increases the 
mean values of Snk(r) and 6Nk{r) after relaxation and the 
temporal fluctuations of both quantities. This supports the ex- 
pectation that with increasing the system size, after relaxation, 
the mean values of Snk (t) and 8N k (r) should become arbi- 
trarily small and so should the temporal fluctuations. 

A different question is what happens to the time scale for 
n(k) to relax to the diagonal ensemble prediction as the sys- 
tem size increases. Previous work on that direction lf30L l3lll 
has suggested a possible intermediate quasisteady regime that 
occurs before full relaxation. This has been seen in recent nu- 
merical work 1I20I1 . and our results for the behavior of n(k) in 
smaller systems may be an indication in that direction. Some- 
thing that is important to keep in mind from our results in this 
section is that different observables may exhibit different re- 
laxation times. In particular, and in contrast to the hardcore 
boson case, the momentum distribution function of fermions 
is an observable that takes a long time to relax to an equi- 
librium distribution. This may have influenced the failure to 



tf in = 1.0, Vf in = 1.0, With 4 

V' in a system with Nf = 7 and L = 21. The initial states were 
selected in the same way as the ones in Fig.[2](T = 3.0). (a) Initial 
and diagonal ensemble results for n(k) and t' = V = 0.32. (b)- 
(e) Time evolution of Sn k for t' = V = 0.32, 0.12, 0.04, and 0.0, 
respectively, (f) Initial and diagonal ensemble results for N(k) and 
t' = V = 0.32. (g)-(j) Time evolution of SN k for t' = V = 0.32, 
0.12, 0.04, and 0.0, respectively. 



observe thermalization in the momentum distribution function 
of a nonintegrable fermionic system in Ref. l23ll . 

It is interesting to note that if V = t' = V = (noninter- 
acting case) the momentum distribution function of fermions 
in a periodic system is a conserved quantity, i.e., for any ini- 
tial state it will not change in time. However, the n(k) of 
hardcore bosons, to which the fermions can be mapped, can 
have a nontrivial dynamics and relaxes to the predictions of 
a generalized Gibbs ensemble, as shown in Refs. |@] and J^]. 
Having V, V' ^ (in our case V = t = 1.0 and V' varies 
between and 0.32) allows the n(k) of fermions to change 
with time (since now the fermions are interacting), but the 
time scale for relaxation to the equilibrium distribution still 
seems to be affected by the fact that in k space the scattering 
between fermions is very special. We find the relaxation time 
for n(k) to be much longer than the corresponding time scale 
for other observables such as N(k). The differences between 
the time scales for the relaxation of different observables (and 
their temporal fluctuations) can be understood in terms of the 
off-diagonal matrix elements of the observables, which will 
be discussed in Sec.lVll 



IV. THERMODYNAMICS 

In Sec. |nil we have shown that, for our systems of inter- 
est, the diagonal ensemble provides an accurate description 
of observables after relaxation. Some observables may take 
longer to relax, and under some circumstances they may not 
even relax at all [see, e.g., n(k) in Fig. |TJe)], but whenever 
we find relaxation to an equilibrium value, it is well described 
by the predictions of Eq. (|9|. This was expected and we re- 
fer the reader to Refs. 114119112511 for similar results in other 
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quantum systems. In what follows, we study how the predic- 
tions of the diagonal ensemble compare to standard statistical 
mechanical ensembles. A great advantage of the infinite time 
average [Eq. is that all time dependence has been removed 
from the time evolution of the quantum-mechanical problem 
and one has a unique result to test statistical mechanics. 

Since the systems on which we have performed the dy- 
namics are isolated, the most appropriate statistical ensemble 
to compare observables after relaxation is the microcanon- 
ical ensemble. As usual, the computations in the micro- 
canonical ensemble are performed averaging over all eigen- 
states (from all momentum sectors) that lie within a window 
[E - AE, E + AE], where E = (^ tm \H fm \ip mi ), and we 
have taken AE = 0.1 in all cases. Similarly to what was done 
in Refs. lfl9ll and i25ll . we have checked that our results are in- 
dependent of the exact value of AE in the neighborhood of 
AS = 0.1. 

The main panel in Fig.|5}a) depicts how the difference be- 
tween the diagonal and the microcanonical ensembles behaves 
as one moves away from integrability. Once again, we find a 
different behavior for n(k) and Nik). For N(k), we find that 
for both effective temperatures [29] considered, the difference 
between the diagonal and microcanonical ensemble is always 
smaller than 1% for t' = V' > 0.06, and one can say that the 
system exhibits thermal behavior. This is very similar to the 
results for the same quantity obtained in the hardcore boson 
systems analyzed in Ref. 12511 . On the other hand, n(k) ex- 
hibits a larger difference between the diagonal and the micro- 
canonical ensemble for the same values of t', V, in particular 
at lower temperatures. For both observables, one can still see 
that as one approaches the integrable point the difference be- 
tween both ensembles increases, signaling the breakdown of 
thermalization in all cases. 

In the inset in Fig. [5} a), we also compare the predictions of 
two diagonal ensembles generated by different initial states, 
which are chosen in such a way that the effective temperature 
of both system after the quench l29ll are the same. The behav- 
ior is qualitatively similar to the one seen in the main panel of 
the same figure and shows that the breakdown of thermaliza- 
tion for both observables is accompanied by a dependence on 
the initial state. 

Results for a smaller system size, with the same density 
and effective temperatures, are shown in the main panel and 
inset in Fig.|5jb). They show that, as expected, the differences 
between the two ensembles for any given value of t' , V' in- 
crease with decreasing system size. Notice that the y axis in 
Figs. [51 a) and|5jb) has a different scale, so that the contrast 
between the results in both panels is stronger than may appear 
at first sight. 

For t' = V' > 0.06, the comparison between the results for 
different systems sizes led us to expect that the predictions of 
the microcanonical ensemble, for both observables, will co- 
incide with the ones of the diagonal ensemble as the system 
size is increased, i.e., that thermalization takes place if one is 
sufficiently far from integrability. We do note that n(k) for 
these fermionic systems is more sensitive to finite-size effects 
than N(k), something that is in contrast to the behavior seen 
for hardcore bosons [25]. 
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FIG. 5: (Color online) (Main panels) Comparison between the pre- 
diction of the microcanonical and diagonal ensembles, for n(k) and 
N(k), as a function of increasing t' , V' for T = 2.0 and T = 3.0. 
Results are shown for (a) L = 24, Nf = 8, and (b) L = 21, N f = 7. 
The diagonal ensembles correspond to the quenches in Figs. Q] [2] 
and [4] (Insets) Comparison between the prediction of diagonal en- 
sembles generated by two different initial states, for T = 2.0 and 
T = 3.0. The results shown are for (a) L = 24, N f = 8, and 
(b) L — 21, Nf — 7. As in the main panels, one of the diagonal 
ensembles is generated by initial states selected from the eigenstates 
of a Hamiltonian with tini = 0.5, Vim = 2.0, the other diagonal 
ensemble is generated by initial states selected from the eigenstates 
of a Hamiltonian with tmi = 2.0, Vini ~ 0.5. The final Hamil- 
tonian (with tfi„ = 1.0, Vfin — 1.0) and the effective temperature 
12911 are identical for both initial states. By relative differences in this 
figure we mean the normalized area between the different ensemble 
predictions for n(k) and N(k). Relative differences are computed in 
exactly the same way as Snk(r) and 5Nk(r) are computed in Eqs. 
d 1 0b and i ll It . respectively. 



For values of t' = V closer to the t' = V = integrable 
point, the outcome of increasing system size is less clear, at 
least for the system sizes that we are able to analyze here. 
Two possibilities emerge, (i) As the system size increases the 
point at which the results for the diagonal an microcanoni- 
cal ensembles start to differ will move toward t' = V' = 0. 
(ii) The very same point will move toward a nonzero value of 
t', V. Here, we are not able to discriminate between those 
two scenarios, and further studies will be needed to address 
that question. 

It is relevant to notice that, for the small systems sizes stud- 
ied in this work, it is important to select the microcanonical 
ensemble when comparing with the outcome of the dynamics. 
This is because finite-size effects can make the predictions 
of different standard statistical mechanical ensembles differ- 
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ent from each other. In Fig. [6] we compare the results of the 
microcanonical ensemble, which we have been using up to 
this point, with the ones obtained with the canonical ensem- 
ble. One can see there that depending on the observable under 
consideration the results of both ensembles can differ up to 
~ 2.5%. The differences are not very much affected by the 
proximity to integrability, and, as expected, they can be seen 
to decrease with increasing system size. One can also see in 
Fig.|6]that the differences between the two ensembles are al- 
ways larger for n(k). Interestingly, they are not as large as the 
ones observed for the same quantity in the case of hardcore 
bosons Il25ll . 
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FIG. 6: (Color online) Relative differences (normalized area) be- 
tween the predictions of the microcanonical and canonical ensem- 
bles for n(k) [upper four curves] and N(k) [lower four curves] in 
systems with L = 21, N f = 7, and L = 24, N f = 8. The rela- 
tive differences are computed in exactly the same way as Sn^ (r) and 
SNk(r) are computed in Eqs. d 1 Ob and i ll It . respectively. 



V. EIGENSTATE THERMALIZATION HYPOTHESIS 

In what follows, we will analyze the reason underlying ther- 
malization in these isolated quantum systems, and the cause 
for the differences seen between n(k) and N(k) when com- 
paring the diagonal and microcanonical ensembles. 

In Sec. IIVI we have shown that away from integrability 
there are regimes in which the system thermalizes, i.e., in 
which the predictions of the diagonal ensemble and the mi- 
crocanonical ensemble coincide. This means that 



Odiag — Omicroci 
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where Me, ae is the number of energy eigenstates with ener- 
gies in the window [E — AE, E + AE]. This result is cer- 
tainly surprising because the left-hand side of Eq. ( TT3T > de- 
pends on the initial conditions through the overlaps of the ini- 
tial state with the eigenstates of the final Hamiltonian (C a ), 
while the right-hand side of Eq. (fT~5b only depends on the en- 
ergy E (as mentioned before, our results do not depend on the 
specific value of AE), 

A possible solution to this paradox was advanced by 
Deutsch l2~lll and Srednicki l22tl in terms of the eigenstate 



thermalization hypothesis (ETH). Within ETH, the difference 
between the eigenstate expectation values of generic few -body 
observables [O aa ] are presumed to be small between eigen- 
states that are close in energy. This implies that if one takes a 
small window of energy AE, as it is done in the microcanoni- 
cal ensemble, all the eigenstate expectation values within that 
window will be very similar to each other. Hence, the mi- 
crocanonical average and a single eigenstate will provide es- 
sentially the same answer, i.e., the eigenstates already exhibit 
thermal behavior. The same will happen with the prediction 
of the diagonal ensemble as long as the |C Q |'s are strongly 
peaked around the energy E. The latter has been shown to 
be the case for quenches in hardcore boson systems in two 
dimensions [19], hardcore bosons in one dimension I25II32I1 . 
softcore bosons in one dimension 13311 . and generic quenches 
in the thermodynamic limit lfl9tl . 




FIG. 7: (Color online) (a),(d) Distribution of |C a | 2 for two of the 
quenches depicted in Fig. [2] for (a) t't in = Vt in = 0.32 and (d) 
t'fin — V'fin = 0.0. In both cases, the final effective temperature 
of the system is T — 3.0 12911 . The vertical dashed lines signal the 
energy E of the time-evolving state. (b),(e) Eigenstate expectation 
values of h(k — 0) in the full spectrum (including all momentum 
sectors) of the Hamiltonian {j} with t = V = 1.0 and, (b) t' = V = 
0.32 and (e) t' = V = 0.0. (c),(f) Eigenstate expectation values of 
N(k — tt) in the full spectrum (including all momentum sectors) of 
the Hamiltonian {]} with t = V = 1.0 and, (c) t' = V' = 0.32 
and (f) t' — V' — 0.0. These results were obtained for systems with 
L = 24 and N f = 8, i.e., the total Hilbert space consists of 735 471 
states. 

In Figs.|7ta) and|7Id), we show the distributions of \C a \ 2 
for two of the quenches studied in Fig. [2] One quench is 
far away from integrability and the other one at integrability. 
There are some features of the distributions of \C a \ 2 that are 
important to mention here, (i) They are neither similar to the 



8 



microcanonical nor to the canonical distributions, (ii) They 
are strongly peaked around the energy E of the time-evolving 
state. Actually, one can see in the figures that the values of 
\C a \ 2 decay exponentially as one moves away from E (sig- 
naled by the vertical line in the figures), (iii) They are not 
related to the effective temperature of the system, which is ex- 
actly the same in both cases (T = 3.0) while the distributions 
of \C a \ 2 are clearly different. The specific values of |C Q | 2 , 
and the exponent of their decay, depend on the initial state 
ll32tl . Because of these properties of the \C a \ 2 distributions, 
they alone cannot explain why thermalization takes place, for 
example, for N(k) when t' = V = 0.32 in Fig.g] 

In the bottom four panels of Fig. [7] we show the eigen- 
state expectation values of h(k = 0) [(b), (e)] and N(k = tt) 
[(c), (f)] for a system far from integrability, with t = V = 1.0 
and t' = V' = 0.32 [(b), (c)], and for a system at integrability, 
with t = V=1.0wdt' = V' = 0.0 [(e),(f)]. Those results, 
together with the \C a \ 2 distributions in Figs. |7J a ) and|7td), 
can help us understand the differences between the diagonal 
and the microcanonical ensembles in Fig. [5ja). 

The most striking feature in the bottom panels in Fig. [7] is 
the contrast between the eigenstate expectation values at in- 
tegrability and far away from it. The former ones exhibit a 
very wide distribution of eigenstate expectation values. This 
is true even if one selects a narrow window at the center of 
the spectrum (equivalent to very high temperatures). As one 
moves away from integrability, one can clearly see that the 
fluctuations between eigenstate expectation values of nearby 
eigenstates reduces dramatically and ETH starts to be valid. 
For t' = V = 0.32, a comparison between the expectation 
values of h(k = 0) and N(k = tt) shows that for the ener- 
gies of the final states created in our quenches, the distribution 
of the eigenstate expectation values of n(k = 0) is relatively 
broader than the one seen for N(k = tt), which explains why 
the differences between the diagonal and microcanonical en- 
sembles in Fig. |3 a) are larger for the former one. It also helps 
in understanding why n(k) is more sensitive to the selection 
of the initial state. 

In order to gain a more quantitative understanding of how 
ETH breaks down as one approaches integrability, we have 
computed the average relative deviation of the eigenstate ex- 
pectation values with respect to the microcanonical predic- 
tion, A mic n(k = 0) and A mic N(k = tt). For any given 
energy of the microcanonical ensemble, these quantities are 
computed as 

A mlc n(k = 0) = E " '"""jf - 0) ~ nm ^ k - 0)1 (16) 



and 



'N(k = tt) = 



n aa( k = 0) 



E„ \N aa {k = tt) - N mlc (k = tt)\ 



E q N aa (k 



(17) 

where n aa (k = 0) and N aa (k = tt) are the eigenstate expec- 
tation values of n(k = 0) and N(k = tt), respectively, and 
n-micik — 0) and N mic {k = tt) are the microcanonical expec- 
tation values at any given energy E. The sum over a contains 



all states with energies in the window [E — AE, E + AE] . As 
discussed before, we have taken AE = 0.1. 

Clearly, for different values of the Hamiltonian parameters, 
the lowest and highest-energy states as well as the level spac- 
ing are different. Hence, in order to make a meaningful com- 
parison between different systems as one moves away from 
the integrable point, we study the behaviors of A mlc n(k = 0) 
and A mlc N(k = tt) as a function of the effective tempera- 
ture T of a canonical ensemble that has the same energy E 
as the microcanonical ensemble. Given an energy E within 
the microcanonical ensemble, the effective temperature can 
be computed by means of Eq. (fT~3T >. We should stress that the 
effective temperature is only used here as an auxiliary tool 
for comparing different systems, i.e., it just provides a unique 
energy scale for comparing observables independently of the 
Hamiltonian parameters, which change the ground-state en- 
ergy and the level spacing. 
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FIG. 8: (Color online) Average relative deviation of eigenstate ex- 
pectation values with respect to the microcanonical prediction vs T 
(see text). Results are presented for: (a),(c) A mlc n(k = 0) and 
(b),(d) A mic N(k = tt) in systems with: (a),(b) L = 24, N f = 8 
and (c),(d) L = 21, Nf = 7. The inset in (d) shows a direct compar- 
ison of A mic N(k = tt) between the systems with L = 24, N f = 8 
(thin lines) and the systems with L = 21, Nf = 7 (thick lines). The 
values of t' , V' in the three cases depicted in the inset are, from top 
to bottom, t' = V = 0, t' = V = 0.06, and t' = V = 0.24, and 
follow the same legends of the main panel shown in (b). In all cases 
t = V= 1.0. 

In Fig. [8] we present results for A mlc n(k = 0) [(a),(c)] 
and A mlc N(k = tt) [(b),(d)] vs T for six different values of 
t' = V' in systems with 24 lattice sites and eight fermions 
[(a),(b)] and 21 lattice sites and seven fermions [(c), (d)]. The 
average relative deviations for both observables are presented 
in the same scale, which immediately allows one to see what 
was already evident in Fig. [7] namely, that the average relative 
deviations for A mlc n(k = 0) are larger (more than two times 
larger) than the corresponding ones for A ml,c N(k = tt) for 
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any given value of t' , V'. For temperatures T > 1.5, one can 
see that those deviations for both observables saturate with in- 
creasing t',V' for t' = V' > 0.12. Comparing the results 
for the same observables but for different systems sizes, one 
can see that with increasing system size the average relative 
deviations are decreasing. This can be better seen in the inset 
in Fig.^d), where we have presented a direct comparison of 
A mlc N(k = 7r) for three different values of i', V' in the sys- 
tems with 21 (thick lines) and 24 (thin lines) sites. At integra- 
bility, one does not see much of a change with changing sys- 
tem size for A mlc N(k — tt), but some reduction can be seen 
for A mlc n{k = 0) [Figs. Eta) andgjc)]. For t' = V > 0, 
we do see a clear reduction in all cases. At low temperatures 
T < 1.5, the density of states is low and in many cases the 
use of the microcanonical ensemble is not justified. From Fig. 
[7] one can see that in that regime (of low energies) the fluc- 
tuations of the observables are very large and thermalization 
is not expected to occur. An important question that needs 
to be addressed in the future is what will happen with that 
window of temperatures (energies) with increasing the sys- 
tem size. For our small systems, we do not see a noticeable 
change with increasing system size. 
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FIG. 9: (Color online) Average relative deviation of eigenstate ex- 
pectation values with respect to the microcanonical prediction at two 
fixed temperatures (T = 2.0 and T — 3.0) and for two system sizes 
(L = 21, N f = 7, and L = 24, N f = 8) vs t', V. Results are 
presented for: (a) A mic n(k = 0) and (b) A mic N(k = tt). In all 
cases t = V = 1.0. 

An alternative way to present some of the results depicted 
in Fig. [8] which can help us make direct contact with the re- 
sults discussed in Fig. [5] is to take two values of the effec- 
tive temperature within the microcanonical ensemble and plot 
the average relative deviation of eigenstate expectation values 
with respect to the microcanonical prediction for those tem- 



peratures as a function of the integrability breaking parame- 
ters t', V. This is done in Fig.|9]for T = 2.0 and T = 3.0, 
and L = 21 and L = 24. Figure [9] clearly shows that the 
breakdown of thermalization seen in Fig.|5]as one approaches 
integrability is directly related to the increase in the relative 
deviation of eigenstate expectation values with respect to the 
microcanonical prediction, i.e., to the increase in the fluctua- 
tions of the eigenstate to eigenstate expectation values of n(k) 
and N(k), or what is the same to the failure of ETH. A com- 
parison between Figs. [9ja) and|9fb) also helps in understand- 
ing why in Fig.|5]the difference between the diagonal and mi- 
crocanonical ensembles is larger for n(k) than for N(k), even 
when one is far from integrability. Figure [9J a) shows that for 
the effective temperatures selected in Fig. [5] A mlc n(k = 0) is 
more than two times larger than A n " c N(k = tt). The reduc- 
tion in A mlc n(k — 0) and A' mlc N{k — tt), with increasing 
system size, can also be seen in Fig. [9] when comparing the 
results for L = 21 and L = 24. 



VI. OFF-DIAGONAL MATRIX ELEMENTS 

While the diagonal elements of the observables, in the ba- 
sis of the eigenstates of the Hamiltonian, helped us understand 
whether after relaxation the expectation values of generic few- 
body observables can be described by standard statistical me- 
chanical approaches, it follows from Eqs. © and (O that the 
off-diagonal elements of the observables can help us under- 
stand the relaxation dynamics and temporal fluctuations of ob- 
servables after relaxation, 

(6(r)) - (OW) = C:C p e^- E ^O a() . (18) 

a/3 

As we pointed out in Sec. Hill one of our findings in this 
work is that the momentum distribution function of fermions 
exhibits a slower relaxation dynamics than the one seen for 
the same observable in an identical system of hardcore bosons 
ll25ll . These differences were seen far from the integrable point 
as well as at the integrable point. In Figs. [Tol a) and [TOl c). 
we compare the off-diagonal elements of h(k — 0) between 
a system of spinless fermions [Fig. [ToT a)l and an identical 
system with hardcore bosons [Fig. [T0l c)l. both systems are 
far from integrability [with t' = V = 0.32] and the central 
eigenstate has an energy for which the temperature of a canon- 
ical ensemble with the same energy would be T — 3.0. It is 
apparent in these two figures that the off-diagonal elements 
of n(k = 0) are larger for the fermions than for the hardcore 
bosons, in particular close to the diagonal. 

The relaxation dynamics of our other observable of inter- 
est N(k), on the other hand, was seen to be very similar 
between the fermions studied in this work and the hardcore 
bosons studied in Ref. i25ll . In both systems, this observ- 
able relaxed very quickly (in a time scale r ~ 1/i). Figures 
[TUl b) and [TUl d) depict results for the off-diagonal elements of 
N(k = tt). The results for fermions and hardcore bosons are 
almost indistinguishable in this case, with all off-diagonal el- 
ements being much smaller than the diagonal ones. The con- 
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FIG. 10: (Color online) Off-diagonal elements of n{k = 0) [(a),(c)] 
and N(k = it) [(b), (d)] for fermions [(a),(b)] and hardcore bosons 
[(c), (d)] far from integrability. Results are shown for the 100 eigen- 
states closest to the one with energy E — —4.26 (T = 3.0, see text) 
for the fermions [(a),(b)] and E = -4.62 (T = 3.0, see text) for the 
hardcore bosons [(c),(d)]. In all cases t = V = 1.0, t' = V = 0.32, 
L = 24, and N f = N b = 8. 




FIG. 11: (Color online) Off-diagonal elements of n(k = 0) [(a),(c)] 
and N(k = ir) [(b), (d)] for fermions [(a),(b)] and hardcore bosons 
[(c), (d)] at integrability. Results are shown for the 100 eigenstates 
closest to the one with energy E = —3.84 (T = 3.0, see text) for 
the fermions [(a),(b)] and E = -3.84 (T = 3.0, see text) for the 
hardcore bosons [(c),(d)]. In all cases t = V = 1.0, i' = V = 0.00, 
L = 24, and Nf = N b = 8. 



trast between the magnitude of the off-diagonal elements of 
h(k = 0) in Fig. [lOta) and N(k = tt) in Fig.GJSb) also ex- 
plains the difference between the relaxation dynamics of both 
observables in fermionic systems depicted in Fig.Q] 

Results for the same quantities for fermions and hardcore 
bosons at the integrable point are shown in Fig. [TT] We find 
the results at integrability in stark contrast with those far away 
from the integrable point. At integrability, the off-diagonal 
elements of both observables are in most cases very small or 
zero, but then there are some states between which the off- 
diagonal elements can be very large (much larger than any 
off-diagonal element seen in the nonintegrable case). In ad- 
dition, very large off-diagonal elements can be seen between 
states that have quite different energies. In Fig.QJJ the largest 
off-diagonal elements in the nonintegrable case always occur 
between states that are close in energies, and they are seen to 
reduce as the energy of the states depart from each other. This 
is not the case at integrability. 

A more quantitative understanding of these issues can be 
gained by taking one state from Figs. [TOl and [TTI and plotting 
the off-diagonal elements between that state (which we take to 
be the one at the center) and all the closest energy eigenstates 
for our two observables of interest. These are depicted in Fig. 
[12] for the nonintegrable case with t' = V — 0.32 [(a)-(c)] 
and at integrability t' = V = [(d)-(f)]. 



First, let us focus on the behavior of n a p(k — 0). Results 
for that quantity in the nonintegrable fermionic and hardcore 
boson cases are shown in panels (a) and (b) for two lattices 
with L = 21 and L = 24, respectively. For each lattice size, 
one can clearly see that the off-diagonal elements of h(k = 0) 
of the fermions are always larger than the ones of the bosons, 
and the differences are largest close to the diagonal. As the en- 
ergies between the states depart from each other n a p{k = 0) 
for bosons and fermions become similar. The comparison be- 
tween the two lattices also shows that with increasing system 
size the off-diagonal elements of the observables decrease, but 
the fact that n a p(k — 0) for the fermions is larger close to the 
diagonal remains. Similarly, we find that with decreasing the 
energy of the state the off-diagonal elements of n a /3(k = 0) 
also become larger, although for temperatures T > 1.0 this 
effect is smaller than the finite-size effect difference between 
panels (a) and (b) in Fig. [12] 

Results for n a p(k — 0) at integrability are presented in 
panels (d) and (e) of Fig. [12] There one can also see that the 
values of n a p(k = 0) are in general larger for fermions than 
for hardcore bosons and that they decrease with increasing 
the system size. However, here fermions and bosons share a 
common feature, namely, large off-diagonal elements can be 
seen between states that have energies that are increasingly 
different from each other. This, in conjunction with the fewer 
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FIG. 12: (Color online) Off-diagonal elements of h(k = 0) [(a),(c)] 
and N(k = ir) [(b), (d)] for fermions [(a),(b)] and hardcore bosons 
[(c), (d)] at integrability. Results are shown for the 400 eigenstates 
closest to the one with energy E = —3.84 (T = 3.0, see text) for 
the fermions [(a),(b)] and E = -3.84 (T = 3.0, see text) for the 
hardcore bosons [(c),(d)]. In all cases t = V = 1.0, t' = V = 0.00, 
L = 24, and N f =N b = 8. 



number of states that have a significant overlap with the ini- 
tial state, can lead to larger temporal fluctuations during the 
time evolution after a quench at integrability. This effect is 
enhanced at lower temperatures where we see some increase 
in the magnitude of the off-diagonal elements and a decrease 
in the number of states that overlap with the initial state (see 
Fig.©. 

Finally, the off-diagonal elements for N(k = ir) far from 
integrability and at integrability are presented in Fig. [12] pan- 
els (c) and (f), respectively. Those panels show that, for the 
density-density structure factor, off-diagonal elements are al- 
ways much smaller than the diagonal ones, and they are rel- 
atively smaller than the ones of ii(k). At integrability [Fig. 
fTZlT )!. one can see a few values of N a p(k = tt) that are 
clearly larger than the rest, but they are still less and relatively 
smaller than the ones seen for n a j3 (k = 0) in Fig. [l"2f e). These 
results explain why N(k) takes less time to relax to the diag- 
onal ensemble prediction and why the temporal fluctuations 
after relaxation are also smaller for this observable. 



VII. SUMMARY 

We have presented a detailed study of the relaxation dy- 
namics after a quantum quench and the description after re- 



laxation of the momentum distribution function n(k) (a one- 
particle observable) and the density-density structure factor 
N{k) (a two-particle observable) of spinless fermions with 
nearest and next-nearest hoppings and interactions in one- 
dimensional periodic lattices. We have also compared some 
results for fermions with those of a similar hardcore boson 
system. We should stress that those models are typical for 
describing the physics of one-dimensional systems. For ex- 
ample, they can be mapped onto the spin- 1/2 XX Z linear 
chain with next-nearest-neighbor interactions. (The observ- 
ables analyzed here are related to the different spin-spin cor- 
relation functions.) Hence, we expect the results reported in 
this manuscript to be generic and apply to other gapless mod- 
els and observables in one-dimensional systems. 

We have shown that, in general, n(k) for fermions exhibits 
a slower relaxation dynamics, i.e., it takes longer to relax to an 
equilibrium distribution than other observables for fermions 
[such as N(k)] and than n(k) for a similar one-dimensional 
system of hardcore bosons (studied in Ref. lf25ll ). Close to and 
at integrability, we have also found that n(k) may even fail to 
relax to an equilibrium distribution while N(k), and n(k) for 
a similar bosonic system, do relax to equilibrium. We have 
shown that this behavior of the dynamics of n(k) is related to 
the particularly large off-diagonal elements of fi(k) in the ba- 
sis of the eigenstates of the Hamiltonian. Those off-diagonal 
elements are larger than the ones for the same observable in 
hardcore boson systems, and much larger than the ones for 
N{k) in the same fermionic system. From the contrast be- 
tween the dynamics of n(k) and N(k) emerges a general pic- 
ture in which some few-body observables in isolated quan- 
tum systems may quickly relax to an equilibrium distribution 
while other observables may take longer to relax, or even fail 
to relax, to an equilibrium distribution. 

We have shown that far from integrability observables af- 
ter relaxation are well described by standard ensembles from 
statistical mechanics, such as the microcanonical ensemble, 
which is particularly relevant to our small and isolated sys- 
tems. The ability of the microcanonical ensemble to predict 
the expectation value of few-body observables after relaxation 
was shown to be related to the validity of the eigenstate ther- 
malization hypothesis, within which the expectation values 
of observables in the eigenstates of the Hamiltonian already 
exhibit thermal behavior. We have also shown that as one 
approaches the integrable point thermalization breaks down, 
with the differences between the microcanonical predictions 
and the results of the relaxation dynamics increasing continu- 
ously as one converges toward the integrable point. We have 
established that there is a clear correlation between the failure 
of the system to thermalize and the failure of the eigenstate 
thermalization hypothesis, i.e., as one approaches integrabil- 
ity the differences between the eigenstate expectation values 
of few-body observables increases between eigenstates that 
are close in energy, and this happens over the full spectrum of 
the Hamiltonian. 

Our results here have been obtained for small one- 
dimensional lattices, which are relevant to current experi- 
ments in one-dimensional geometries B3411 . However, several 
important questions remain to be addressed in future works. 
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For example, what happens to the point at which thermaliza- 
tion breaks down (at all temperatures) as one increases the 
system size. Far from integrability, we find that for our small 
systems thermalization occurs whenever the effective temper- 
ature is T > 1.5 M29I1 . Another question that needs to be ad- 
dressed is what happens with the window of energies where 
thermalization fails to occur far from integrability as one in- 
creases the system size. Finally, there is the question of the 
time scale that n(k) of the fermions takes to relax to the ther- 
mal distribution, and the emergence of intermediate quasista- 
tionary distributions lf3QL l3ltl . as the size of the system is in- 
creased. Experiments with ultracold gases will generate many 



new questions and help address the current ones. 
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